Giant anisotropy of Zeeman splitting of quantum confined acceptors in Si/Ge 
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Shallow acceptor levels in Si/Ge/Si quantum well heterostructures are characterized by resonant 
tunneling spectroscopy in the presence of high magnetic fields. In a perpendicular magnetic field 
we observe a linear Zeeman splitting of the acceptor levels. In an in-plane field, on the other hand, 
the Zeeman splitting is strongly suppressed. This anisotropic Zeeman splitting is shown to be a 
consequence of the huge light hole-heavy hole splitting caused by a large biaxial strain and a strong 
quantum confinement in the Ge quantum well. 

PACS numbers: 71.18.-|-y, 71.70.Fk, 73.21.Fg 



Spintronic and quantum computing are novel de- 
vice concepts relying on quantum mechanical coherence. 
Si/Ge-based systems are promising candidates offering 
long spin coherence times fast operations, and a 

well-established record of scalable integration. These im- 
portant properties are also crucial requirements |^ IH 1^ 
for implementing multi-qubit operations in a future quan- 
tum computer. One concept that may form the techno- 
logical basis of a quantum computer is the spin-resonance 
transistor (SRT) 0. Vrijen et. al. ^ proposed an SRT 
where the electron spin manipulation is realized using the 
change in gr-factor between Si-rich and Ge-rich environ- 
ments of a Si/Ge- heterostructure. However, engineering 
the g-factor in such systems is complicated by the fact 
that the electron states in Si are in the X valleys whereas 
in Ge the electrons are located in the L valleys . This 
problem does not arise for the valence band states, as 
both Si and Ge have their valence band maximum at 
the r point. Thus, valence band states in Si/Ge are a 
promising choice for g-factor engineering in search for 
spin manipulation. 

In this paper we have analyzed the (7-factor of shallow 
acceptor levels in a Si/Ge-heterostructure by resonant- 
tunneling spectroscopy. We find that their effective 
(jf-factor is highly anisotropic, giving a large Zeeman 
splitting of the acceptor states in a perpendicular field, 
whereas we cannot resolve any Zeeman splitting in in- 
plane fields up to 18 T. This giant anisotropy provides 
the possibility to tune the coupling of the holes to an 
external magnetic field by a gate-controlled shift of the 
hole wave function |loj in spintronic devices. 

For a proper understanding of acceptor levels it is es- 
sential to take into account the fourfold degeneracy of 
the valence band at the F point (Fig. ^ which reflects 
the fact that the bulk valence band edge in these mate- 
rials is characterized by an effective angular momentum 
J — 3/2 [ill H^ . As the symmetry of the crystal is 
reduced due to biaxial strain and a confinement poten- 
tial, the degenerate states split into heavy-hole (HH) sub- 
bands with J2 = ±3/2 and light-hole (LH) subbands with 
Jz = ±1/2. Here, the quantization axis for the angular 
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FIG. 1: Qualitative sketch of the valence band states at differ- 
ent perturbation conditions. Already at B = 0, biaxial strain 
and a confinement potential reduce the fourfold degeneracy to 
twofold, giving heavy-hole (HH) states with Jz = ±3/2 and 
light-hole (LH) states with Jz — ±1/2. Depending on the 
orientation of the magnetic field, the twofold degeneracies of 
the HH and LH states will be lifted or preserved. 



momentum is the z-axis perpendicular to the epitaxial 
layer. So both parameters, the confinement potential 
and the built-in strain, substantially influence the energy 
levels of an acceptor in a quantum well (QW) [isl Il4l |. 
In a magnetic field B± orientated perpendicular to the 
epitaxial layer we get a Zeeman splitting of HH and LH 



states. 



^^HH (LH) 



5hh (lh)/^bBx, where g^^ (^h) is the 



g-factor of the HH (LH) states in a perpendicular field 
and /iB is the Bohr magneton. But for an in-plane mag- 
netic field -By the linear Zeeman splitting of HH states 
is suppressed because there is no By-induced direct cou- 
pling between these states, (HH|J • By|HH) = 0, where 
J is the vector of J = 3/2 spin matrices HI, E lUlIll ■ 
This does not apply for LH states which show a signif- 
icant Zeman splitting AiJ^jj > 0. We emphasize that 
the vanishing Zeeman splitting of HH states in a paral- 
lel field reflects the fact that the HH-LH splitting in our 
samples is much larger than the maximal Zeeman ener- 
gies (- 7 meV) 0. 

Figure[21shows the layer sequence and the valence band 
profile of our samples. They are prepared by growing a 
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FIG. 2: (a) Layer sequence of the heterostructure. (b) Self- 
consistently calculated valence band profile of the active re- 
gion. The solid line shows the shape of the heavy-hole (HH) 
subband and the dashed line the light-hole (LH) subband. 
Due to biaxial strain and a strong quantum confinement a 
huge HH-LH splitting results. 



100 nm thick Si:B layer (p = 8 x 10^^ cni-^) on a (DDI) 
p+-Si substrate. The active region of the samples consists 
of two 10 nm thick SiGe QWs separated by a 10 nm thick 
Si barrier. In the center of the Si barrier a 4 monolay- 
ers (ML) thick Ge QW is embedded. Finally, the active 
region is capped with 100 nm Si:B {p — 8 x 10^^ cm~'^). 
For the dc-transport measurements we have fabricated 
diodes with lateral diameters of 1 /xm. Measurements 
were performed at temperatures down to T = 50 mK 
and using magnetic fields up to 18 T. 

Figure Ola) shows a typical current-voltage (/ — V) 
characteristic of a diode at T = 50 mK. A staircase struc- 
ture is observed that is even better resolved in the dif- 
ferential conductance (dl/dV) curve shown in Fig. |3Ib). 
The steplike increase of the current is in a bias range 
which is about 300 mV below the onset of resonant tun- 
neling of holes through the 2D subbands of the central Ge 
QW so that this mechanism cannot explain the current 
steps. We attribute these current steps in the pA-range 
to tunneling processes of holes through zero-dimensional 
acceptor levels of Boron dopant-atoms which have mi- 
grated into the Ge QW from the highly doped Si:B con- 
tact regions. 

A resonant tunneling process through an acceptor level 
Eg occurs each time Eg is in resonance with the Fermi 
energy Ep of the emitter. The bias position of a cur- 
rent step is given by Vs — (Eg — Ep)/ae, where a is the 
bias-to-energy conversion coefficient. We determine a 
from the temperature-dependent broadening of the cur- 
rent step edges. As a measure of this broadening, we 
use the full width at half maximum of the corresponding 
differential conductance peaks, AFg [see Fig. Efb)]. It 
increases according to (aeAI^)^ = [l^Esf + (3.53 fer)^ 
where the term 3.53 fcT stems from the broadening of 
the Fermi function characterizing the carrier distribution. 
Using this equation we obtain a = 0.5 ± 0.1 for several 
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FIG. 3: (a) Current-voltage (/ — V) characteristics of a diode 
with a diameter of 1 ^m. at T = 50 mK and (b) the corre- 
sponding differential conductance (dl/dV). A current step 
respectively a differential conductance peak occurs whenever 
resonant tunneling through a shallow acceptor level in the Ge 
quantum well is energetically possible. 




FIG. 4: Gray-scale plot of the differential conductance at 
50 mK for a magnetic field orientated perpendicular to the 
epitaxial layer, where a dark (bright) shade corresponds to 
small (large) conductance. Exemplarily, the top graph shows 
the splitting of a level, marked by an arrow at 6.5 mV. 



peaks in both bias polarities. 

Figure2]shows a gray-scale plot of the differential con- 
ductance dl/dV as a function of an external magnetic 
field orientated perpendicular to the epitaxial layer, B^. 
In the voltage range from —15 mV to 15 mV the conduc- 
tance maxima exhibit a linear splitting as a function of 
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FIG. 5: Gray-scale plot of the differential conductance at 1 K 
for a magnetic field orientated parallel to the epitaxial layer, 
where a dark (bright) shade corresponds to small (large) con- 
ductance. The right arrow points to the same level as in 
Fig.EI 



B±. All levels show accurately the same splitting, as in- 
dicated by the parallel evolution of the conductance max- 
ima. As an example, the upper panel of Fig. 01 shows the 
splitting of the conductance maximum at 6.5 mV. The 
gradient of the splitting is dAVz/dS^ = 0.73 mV/T. We 
attribute the linear splitting to the Zeeman splitting of 
the HH sublevels, AEz = g^YifJ,BB±. The (/-factor 
can be determined using 



ffHH 



ae dAVz 
Mb dS^ 



(1) 



With a = 0.5 and dAFz/dB^ = 0.73 mV/T we obtain 
f/j^jj = 6.3. This value agrees well with optically mea- 
sured 5-factors |0, of group-HI impurities such as B 
in Ge. This confirms the assumption that the observed 
levels belong to Boron dopant-atoms which have diffused 
from the heavily doped contact regions into the Ge QW. 

Using the B± dependence of the conductance we can 
obtain an upper limit for the radial extent A of the 
wave function of a hole bound to an acceptor level. 
According to first-order perturbation theory, a ground 
state acceptor level is affected by a diamagnetic shift 
A£;d ~ e^BlX^/8m*, with m* « 0.28mo the HH ef- 
fective mass of Ge. (We neglect here the modifications 
of the hole wave functions caused by the confinement in 
the narrow QW.) But even at B± = 18 T a diamag- 
netic shift of the levels cannot be observed in our mea- 
surements, only the contribution of the linear Zeeman 
splitting can be seen in Fig.^] This implies that the dia- 
magnetic shift Ai?D of the acceptor level is smaller than 
the width AEg ~ 0.15 meV of the conductance peak at 
50 mK so that A must be smaller than 2.5 nm. There- 
fore, A of the acceptor wave functions in the thin Ge QW 
is of the same order as the QW width. We remark that 
Bastard derived a simple model to estimate A for an 



impurity in a narrow QW which yields for our system pa- 
rameters A « O.Sflg, where Og ~ 3.1 nm is the effective 
Bohr radius for heavy holes in Ge. 

Next we present in Fig. (Sjour results for the measured 
conductance in an in-plane magnetic field . While we 
saw in Fig. ^ that B± gives rise to a significant Zeeman 
splitting of the acceptor levels linear in B±, it is most 
remarkable that up to 18 T most conductance maxima 
are not at all influenced by an in-plane magnetic field -B|| . 
The conductance maximum at 6.5 mV, which exhibits a 
pronounced linear splitting for B± (arrow in Fig.^, does 
not show any splitting in the case of an in-plane magnetic 
field (arrow in Fig. [5)1. If a splitting exists, it must be 
smaller than the width of the conductance peak which is 
about 0.35 meV 22] here. 

The giant anisotropy of the Zeeman splitting is a con- 
sequence of the effective spin J — 3/2 of the valence 
band states (Fig.QJ. For a detailed interpretation of our 
experimental results, we have performed self-consistent 
calculations in the multiband envelope-function approxi- 
mation of the valence band profile of the active region 
using the nominal growth parameters. The results of the 
calculation are plotted in Fig. I^^b). The solid (dashed) 
lines show the HH (LH) subbands and the strain-split 
effective potentials for these states. The calculation pre- 
dicts a splitting of the lowest HH and LH subbands of 
about 200 meV. This huge HH-LH splitting is caused 
by the strong quantum confinement of the thin Ge layer 
and the large biaxial strain due to the lattice mismatch 
between Ge and Si. 

The behavior of the HH states in our device is in sharp 
contrast to electron states for which it is well-known that 
the Zeeman splitting is proportional to the total magnetic 
field B irrespective of the orientation of B relative to the 
epitaxial layer. Furthermore, confinement potential and 
strain do not affect the Zeeman energy of electron states. 
In the case of HH states, on the other hand, the Zeeman 
splitting in a field By competes with HH-LH splitting; 
Zeeman splitting is the smaller the larger the HH-LH 
splitting and vice versa The appropriate situation 
can be created in a narrow QW or by application of uni- 
axial or biaxial stress. Our samples satisfy both of these 
requirements so that we obtain a huge HH-LH splitting, 
as can be seen in the band profile in Fig. Elb), result- 
ing in a vanishingly small Zeeman splitting in a field -B|| . 
This explains why we do not observe a Zeeman splitting 
in Fig. 

In the discussion of Fig.^only the isotro pic p art of the 
bulk Zeeman term was taken into account |ll[ Il2l | . This 
is the dominant term. The anisotropic part is typically 
two orders of magnitude smaller than the isotropic part 
and the calculations predict for our structure that it gives 
rise to a linear splitting with AS^jj w 0.18 meV at B\\ = 
18 T. Such a small splitting cannot be resolved in our 
experiment due to the width of the conductance peaks. 



4 



It corresponds to g'^^ « 0.17 which is almost two orders 
of magnitude smaller than g^^i — 6-3 

In rare cases of conductance peaks we see a slightly 
different behavior. In Fig.Eltwo conductance maxima at 
+10 mV and —10 mV indicate a small nonlinear splitting 
above 10 T. If HH-LH coupling is taken into account, we 
get a splitting cubic in B\\ and inversely proportional to 
the HH-LH splitting, A£;z oc B\l\EYi'i^~EiM\ (Ref. [l^). 
For a fully strained system [Fig.|2tb)] the calculated split- 
ting due to this term is even smaller than the splitting 
due to the anisotropic Zeeman term. However, it is con- 
ceivable that the levels showing a splitting nonlinear in 
are related to shallow acceptors situated in sample re- 
gions of slightly relaxed strain (e.g., close to misfit dislo- 
cations) . In these regions the HH-LH splitting is thus re- 
duced and the cubic Zeeman splitting increases for these 
levels. This can explain why a nonlinear splitting is ob- 
servable for the two conductance peaks at ±10 mV, but 
not for the majority of the conductance resonances which 
are due to impurities in highly strained regions. 

In conclusion, we have performed a detailed study of 
Zeeman splitting of shallow acceptor levels in a thin 
Si/Ge/Si quantum well, by using resonant tunneling 
spectroscopy. In a magnetic field orientated perpendic- 
ular to the layer a large linear Zeeman splitting can be 
observed for magnetic fields up to 18 T. In an in-plane 
magnetic field the Zeeman splitting is suppressed. The 
giant anisotropy of the Zeeman splitting is a consequence 
of the huge heavy hole-light hole splitting produced by 
a large biaxial strain and a strong quantum confinement 
in the narrow Ge quantum well. It opens a new way to 
(^-factor engineering for spintronics and quantum com- 
puting. 



[1] 



Present adress: Department of Physics, Northern Illinois 
University, De Kalb, IL 60115 

D. D. Awschalom, D. Loss, and N. Samarth, Semiconduc- 
tor Spintronics and Quantum Computations (Springer, 
Berlin, 2002). 



[2] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 
(1998). 

[3] A. M. Tyryshkin, S. A. Lyon, W. Jantsch, and 

F. Schaffler, Phys. Rev. Lett. 94, 126802 (2005). 
[4] A. M. Tyryshkin, S. A. Lyon, A. V. Astashkin, and A. M. 

Raitsimring, Phys. Rev. B 68, 193207 (2003). 
[5] J. Preskill, Proc. R. Soc. Lond. A 454, 385 (1998). 
[6] R. de Sousa and S. Das Sarma, Phys. Rev. B 67, 033301 

(2003); R. de Sousa and S. Das Sarma, Phys. Rev. B 68, 

115322 (2003). 

[7] B. E. Kane, Nature 393, 133 (1998); A. J. Skinner, 
M. E. Davenport, and B. E. Kane, Phys. Rev. Lett. 90, 
087901 (2003). 

[8] R. Vrijen, E. Yablonovitch, K. Wang, H. W. Jiang, A. Ba- 
landin, V. Roychowdhury, T. Mor, and D. DiVincenzo, 
Phys. Rev A 62, 012306 (2000). 
[9] F. A. Baron, A. A. Kiselev, H. D. Robinson, K. W. 
Kim, K. L. Wang, and E. Yablonovitch, Phys. Rev. B 
68, 195306 (2003). 

[10] I. Malajovich, J. J. Berry, N. Samarth, and D. D. 
Awschalom, Nature 411, 770 (2001). 

[11] J. M. Luttinger, Phys. Rev. 102, 1030 (1956). 

[12] R. Winkler, Spin-Orbit Coupling Effects in Two- 
Dimensional Electron and Hole Systems (Springer, 
Berlin 2003). 

[13] W. T. Masselink, Y.-C. Chang, and H. Morkoc, Phys. 

Rev. B 32, 5190 (1985). 
[14] V. Ya. Aleshkin, V. I. Gavrilenko, D. B. Veksler, and 

L. Reggiani, Phys. Rev. B 66, 155336 (2002). 
[15] H. W. van Kesteren, E. C. Cosman, W. A. J. A. van der 

Poel, and C. T. Foxon, Phys. Rev. B 41, 5283 (1989). 
[16] E. Glaser, J. M. Trombetta, T. A. Kennedy, S. M. Prokes, 

O. J. Glembocki, K. L. Wang, and C. H. Chern, Phys. 

Rev. Lett. 65, 1247 (1990). 
[17] S. Y. Lin, H. P. Wei, D. C. Tsui, J. F. Klem, and S. J. 

Allen, Jr., Phys. Rev. B 43, 12110 (1991). 
[18] G. Goldoni and A. Fasolino, Phys. Rev. B 48, 4948 

(1993). 

[19] H. P. Soepangkat and P. Fisher, Phys. Rev. B 8, 870 
(1973). 

[20] J. Broeckx, P. Clauws, K. van den Steen, and J. Vennik, 
J. Phys. C: Sohd State Phys. 12, 4061 (1979). 

[21] G. Bastard, Phys. Rev. B 24, 4714 (1981). 

[22] The transport measurements in Fig. |S]were performed at 
1 K which resulted in a full width at half maximum of 
the corresponding conductance peaks of 0.35 meV. 



